Electron-Dephasing Time in A Two-Dimensional Spin-Polarized System with Rashba 
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We calculate the dephasing time i>(B) of an electron in a two-dimensional system with a Rashba 
spin-orbit interaction, spin-polarized by an arbitrarily large magnetic field parallel to the layer. 
t$(B) is estimated from the logarithmic corrections to the conductivity within a perturbative ap- 
proach that assumes weak, isotropic disorder scattering. For any value of the magnetic field, the 
dephasing rate changes with respect to its unpolarized-state value by a universal function whose 
parameter is 2Ez / Esoi {Ez is the Zeeman energy, while Esoi is the spin-orbit interaction), con- 
firming the experimental report published in Phys. Rev. Lett. 94, 186805 (2005) . In the high 
field limit, when 2Ez » Esoi, the dephasing rate saturates and reaches asymptotically to a value 
equal to half the spin-relaxation rate. 
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PACS numbers: 

INTRODUCTION 

In two-dimensional electron gases (2DEG) with spin- 
orbit interaction (SOI), the skew-scattering of the elec- 
tron spin leaves a trademark fingerprint on the quantum 
corrections to the conductivity in the form of a positive 
contribution, the antilocalization term. This results from 
the SOI mediated coupling of trajectories of electrons 
with opposite momenta and opposite spins, a configura- 
tion aptly named a singlet Cooperon, in the impurity- 
averaged diffusion equation [l[. The phase difference in 
the interfering paths of the two electrons is measured by 
the dephasing time 

A magnetic field of intensity B applied parallel to the 
2D layer diminishes the spin-scattering effect of SOI by 
aligning the spins parallel to its direction. In this situ- 
ation, t^(B) is expected to depend on two relevant pa- 
rameters: the Zeeman splitting, proportional with the 
effective gyromagnetic factor 7, Ez = 2^B and the spin- 
orbit energy, expressed in terms of the spin-relaxation 
timer s (0) Q, E SO i = h/r s (0). 

A theory of this effect, developed at low magnetic fields 
by considering the Zeeman interaction as a perturbation 
on the electron diffusion equation [H, finds that, for a 
SOI with linear (Rashba [H) momentum coupling, t^B) 
is proportional to {Ez / Esoi) 2 [4]. This behavior has 
been confirmed by several separate experiments 0, Q . 

The linear {Ez / Esoi) 2 dependence of i>(B) inferred 
from the low field estimates has been challenged recently 
by data obtained by Meijer et al. Q for intense magnetic 
fields. When Ez/Esoi >> 1 the dephasing time satu- 
rates, becomes independent of such system parameters as 
the electron density, Rashba splitting, or the elastic scat- 
tering time, and is represented by a universal function of 
{Ez / Esoi) 2 ■ 

Challenged by these experimental findings, we calcu- 
late the dephasing time and study its variation with Ez 
from a general theory of the localization effects in a spin- 



polarized 2DEG with Rashba spin-orbit coupling. The 
first important result of this paper is that within the 
weak, isotropic scattering approximation, r^(-B) is, in- 
deed, a universal function of {2Ez / Esoi) 2 ■ This de- 
pendence, rather than the {Ez / Esoi) 2 parametrization 
discussed in the experiment, might be explained by the 
argument that the additional dephasing introduced by 
the longitudinal magnetic field in the trajectories of the 
opposite spin electrons that form the singlet Cooperon is 
such that it preserves the relative orientation of the two 
spins, making 2Ez (corresponding to two spin flips) the 
appropriate energy scale. The second important result is 
that in the high-field limit, the calculated relative change 
of the dephasing rate in respect to its unpolarized value 
tends asymptotically to 1/2t s (0). Of course, at low fields 
we recover the dephasing rate obtained in Ref. |4J. 

The object of this study is a 2DEG, spin-polarized by 
a magnetic field parallel to the layer, B = Bz. The mag- 
netic field determines an imbalance between the num- 
ber of spins parallel to the field, n\ and those oppo- 
site to the field, nj.. The degree of spin polarization, 
£ = (n-f — ni)/{n^ + 71 jj, varies continuously from — 1 
to 1 as a function of B and is considered a parameter 
of the problem. The direction of the magnetic field, in 
the plane of the gas, corresponds to the quantization axis 
of the electron spin, z {x is in the plane, while y is per- 
pendicular on the plane). The non-interacting electron 
energies, in the absence of the spin-orbit interaction, re- 
flect the Zeeman splitting and are written, for an electron 
of momentum p, spin er, and effective mass m* (consid- 
ered to be spin independent), as eg >a = p 2 /2m* — jcrB, 
where a = 1 for a spin parallel to the field and a = — 1 
for the opposite. 

In the presence of the Rashba spin-orbit coupling, the 
single-electron Hamiltonian is 



Hsoi — 



p 

2m 



X . 



(fxy)-ja-B 



(1) 



2 



with A the coupling constant of the Rashba interaction, 
and a the Pauli spin operator. 

To calculate the dephasing time, we generalize the well 
known theory of localization Q for a spin-polarized sys- 
tem. By considering the spin polarization as a preexist- 
ing condition of the problem, rather than a perturbation 
in comparison with SOI, as in Ref. we allow an un- 
constrained relationship between Ez and Esoi to occur. 
The logarithmic corrections to the conductivity are gen- 
erated by the poles of a quantum diffusion equation, or, 
in an equivalent representation [9(, by the eigenvalues of 
the Cooperon equation. The propagator (or Cooperon) 
represents the convolution in the momentum space of two 
single-particle Green functions, averaged over the impu- 
rity configuration. To construct the propagator, we start 
by first obtaining an expression for the single-electron 
Green function in a spin-polarized electron system with 
Rashba interaction. 

In the absence of SOI, the single particle Green func- 
tion in the 2DEG is represented by a diagonal 2x2 matrix 
in spin space: 



G R ' A (p,u) = 



n R,A 

(jr.* 



(p>) 









Gf' A (p,v) 



(2) 



where G^' A (p,uj) = (lu — TC ± -^r)^ 1 ar e the spin- 
dependent retarded (R) (with +) and advanced (A) (with 
— ) components. The impurity scattering is assumed to 
be isotropic, characterized by a rate tq = 2ttNqu 2 , with 
N = m/2irh 2 the single spin density of states at the 
Fermi surface, and u 2 the mean-square impurity poten- 
tial, considered spin- independent. 

The non-averaged propagator, P°(q, SI) defined by 

P°(q, Sl) = [duJ2 ")G A (-p + q, u + Sl), (3) 

p 

involves two Green functions associated with two differ- 
ent spin 1/2 particles. Its matrix representation occurs in 
the Si ® S2 space, where the Pauli spin operators and the 
identity that act as a basis are, respectively, {Ii,<J x ,y,z} 
and {I2, Px,y,z}- In this description, P°(q,Sl) becomes: 



P°(q,Sl) = ]T 



+P° (q,Sl) 



h 



2 

crp~z 



■(4) 



The coefficients P° a , (q, SI) are calculated from Eq. ([3]) 
for the corresponding pairs of spins in the perturbative 
approach of the weak scattering approximation, by con- 
sidering 1/to 3> much larger than all the other frequen- 
cies involved: 



2tt7V t [1 
2ttN t [1 



iSlr - D a q 2 T ] , (5) 
• i(Sl + 2a~fB/h)T - D a q 2 T } (6) 



D a = Dq(1 + er£) are the spin dependent diffusion coef- 
ficients, which in the case of the unpolarized gas become 
equal to the diffusion coefficient in 2D, D = v f Tq/2 , 
where vp is the Fermi energy. 

The introduction of SOI changes substantially this pic- 
ture because the spin direction is no longer parallel to 
the applied magnetic field, but rather undergoes a pro- 
cess of randomization induced by the coupling to the or- 
bital motion. A single-electron Green function that in- 
corporates both interactions can be obtained by using 
the exact eigenfunctions and the corresponding eigenval- 
ues of Hsoi- Because of the spin mixing, however, this 
algorithm is laden with serious mathematical difficulties, 
especially where the impurity average is concerned. 

In a different approach, adopted here, we exploit the 
linear coupling of the electron momentum to the spin in 
the Rashba term. This coupling allows the introduction 
of a spin-dependent vector potential m\{u xy)/h that is 
being used as a generator of a non-abelian unitary trans- 
formation in the spin space to produce the Green's func- 
tion of the spin-polarized electrons in the presence of the 
spin-orbit interaction. Such an approximation has been 
discussed before in connection with a magnetic vector 
potential ||, with the Aharonov-Casher effect and 
with SOI in parallel quantum dots Consequently, 
we define 



(7) 



as the single particle Green function in the presence of 
SOI, with A = Xm/h a reduced Rashba coupling con- 
stant. On account of SOI, the new Green function Q 
is not diagonal. Accordingly, the non-averaged propaga- 
tor V°(q, SI), defined by Eq. ([3]), becomes a sixteen term 
sum in the vector space Si (g> S2, that reflects the skew- 
scattering of the electron spins. 

Performing the statistical average over the impurity 
configuration requires the propagator V(q, St) to satisfy: 



V{q, St)=u 2 + u 2 P°(q, Sl)V(q, SI) 



(8) 



An elegant solution to Eq. ([5]) is obtained in a basis of 
the 5*i ® 5*2 vector space in which the the kernel P°(q, SI) 
assumes a diagonal form. The vectors Vt'i(r) of this basis 
are found by solving 



dr'V°(f,r')^i(r') = A 4 * 4 (r) 



(9) 



for the corresponding eigenvalues A$. In this representa- 
tion, P°(f,P) = J2i A-i{q,Si)^i{f)^(P) leading to an 
immediate solution to Eq. ([5]): 



u 2 A t {q,Sl) 



^•(r)**(r') . (10) 



The logarithmic corrections to the conductivity are then 
given by the poles of P(q,Sl), estimated in the limit of 



3 



small q: lim^o(l — u 2 Ai) [8(. The problem is centered, 
therefore, on solving the eigenfunction-eigenvalue equa- 
tion of the non-averaged propagator. 

Eq. (J9]) is linearized by expanding P® a , and \&i(r') 
up to second order around r Q. The final result in- 
volves the Fourier transforms of the propagators in the 
absence of SOI, P® a , , defined in Eq. ([6]), and their deriva- 
tives with respect to q. Simultaneously, the eigenfunc- 
tion %(P) is replaced by = ^i(r) + V* i (r)Af + 

iV 2 * J (r)(Ar) 2 . *i(f) is a spinor in Si ® S 2 whose gen- 
eral form is a linear combination of eigenvectors of S Z i, 
in Si ® S2, \S, >i for spin up and |ry >j for spin down 
(t = 1,2) : *(r) = e?"(a|£ >i |£ >2 >i \v >2 

+c|?7 >i |£ >2 +d|?7 >i |?7 >2) , where a,b,c,d are nu- 
merical coefficients. With this test solution, the charac- 
teristic equation for the eigenvalues of the propagator is 
obtained to be 



Zi] q z q z -A 

q z Z | —A q z 

q z -A Zoi q z 

-A q z q z Z n 



0, 



(11) 



where the diagonal coefficients that in 
corporate the eigenvalues A are Zi a 
[l + iVLt - D a TQ (q 2 - AXq x - 6A 2 ) - A] /2XD a T 
and Zo a = 

[l+i{Sl + 2o 1 B/K) T - L> ff to (q 2 - 2A 2 ) - A] /2AD ff To . 
q x and q z are the corresponding components of q. 

The straightforward evaluation of the determinant 
leads to a quartic equation in A (imbedded in Zi a ): 

(Zo^Zqi - A 2 ) (Zi-j-Zix - A 2 ) = 

q 2 z (Z 11 + Z n + 2A) (Z 0T + Z Ql + 2A) . (12) 

Even though, in principle, this equation can be solved 
for any values of q, for the problem at hand suffices to 
obtain its solutions for q = 0. It is also instructive to 
investigate the solutions of this equation in the absence of 
the magnetic field. In the relevant limit q — ► 0, we obtain 
Z\± = ±A and Zq± — ±A, which generate the following 
values for the logarithmic corrections to the conductivity: 



lim ■ 



TO 



-iil + 8A 2 A) 

-in + 4A 2 A) 

-in + 4A 2 L>o 

-in 



(13) 



By using the definitions of the diffusion coefficient and 
of the reduced Rashba interaction, we identify ADqX 2 = 
2X 2 kpT (kp is the Fermi momentum) as the Dyakonov 
spin relaxation rate, 1/t s (0) [2(. This is associated with 
the SOI induced rate-of-change of the in-plane spin com- 
ponents : dSi/dt = —Si/ru, (i = x, z), where tu = t s (0). 
The rate of change of the perpendicular spin component, 
dSy/dt is described by r yy = t(0)/2. The imaginary fre- 
quency is replaced by the dephasing rate, r7 (0). With 



these substitutions, we readily regain the traditional form 
of the conductivity poles expressed in terms of the spin- 
relaxation times as derived in Ref. [l[ . One can also show 
that the eigenfunction fy(r) associated with the eigen- 
value tT^O) corresponds to the pairing of two electrons 
of opposite spins and momenta, i.e. the singlet Cooperon, 
and generates the antilocalization correction [l[. 

Inspired by this analysis, we rewrite Eq. (|12p in terms 
of a new set of unknowns, = Zi^ a + A, such that 



"IT + U H 



- A 



u 01 u 0j 



- A 



(14) 



In this form, the limit q — > can be taken and the re- 
sulting two uncoupled quadratic equations can be solved 
independently. With \i\ a = 27riVoTo[l + ifiro — 4A 2 £) (T t ] 
and Li Qcr = 2ttN [1 + i(n + 2jaB/h)T - 2X 2 D t ], the 
solutions are: 



A 



1 r 

2 [ 



f i ll -2X 2 (D^+D l )T ± 



(15) 



l( M - Mi;) - 2A 2 t 2 (D t - D{)] 2 + \QX^D^D i 



= A quick inspection shows that 1 — u 2 Ai t ± \^ =0 /tq and 
1 — w 2 Ao,- |^=o / T o correspond to the first three lines of 
Eq. (113|) . Their dependence on the magnetic field is pro- 
ducing a negligible effect on the logarithmic corrections 
to the conductivity compared with the unpolarized case. 

The opposite is true, however, about the effect of the 
magnetic field on the solution that describes the weak 
antilocalization, 



lim 

q^0 



1 - m 2 A ,- 

T0 



-ifi + 2D A 2 



2A)A 2 W1 + 2C 



vyB 



KD t X 2 



■yB 



hD Q X 2 T Q 



.(16) 



In analogy with the unpolarized case, described by the 
last line of Eq. (fl3|) . a dephasing time t^B) is defined as 
a measure of the antilocalization correction in the pres- 
ence of a magnetic field: 



2DnX 2 




( ijB V 
U-DoroAV 



(17) 



[5fte designates the real part of the expression.] 

Following the notations of Ref. 0, we introduce the 
variation of the dephasing rate from the unpolarized 
case, T S AB) = [t^B)]^ 1 — [t^(O)]^ 1 and recognize that 

(jB/hDoroX 2 ) 2 = {4 7 B/[?i/t s (0)]} 2 = (2E z /E SO i) 2 . 
Moreover, since on account of the large electron density 



4 



and the low effective mass the spin polarization is very 
small, in the following considerations £ is set to 0. Thus, 
by extracting the real part of Eq. (fTT|) we obtain that the 
change of the dephasing time induced by the longitudinal 
field is simply: 

r.(o)r;(B) = 



1 

- < 1 



\ 











1 - 


( 2 M 2 + 


1 - 


(2 Ez y 


\Esoi J ^ 


\Esoi J 



(18) 



We compare the relative dephasing rate give by Eq. (fl8|) 

with tl 1 ^ ovnarimDnf /lofa in TT'irr / j 1 [ \ 




O O o o%> 



° % 



o experimental data 
— theoretical fit 



I 
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FIG. 1: The calculated dephasing time from Eq. 
pared with the experimental data of Meijer et al. 



Tap is com- 
for C = 0. 



The dephasing mechanism suggested by Eq. (fTg)) can 
be understood by considering that in the presence of the 
longitudinal magnetic field the down-spin component of 
the singlet Cooperon is forced to undergo a spin-flip. 
This action is counterbalanced by the spin-orbit coupling 
which will redirect the spin in some other direction. For 
the singlet Cooperon to continue to exist and thus con- 
tribute to the antilocalization, the dephasing of the orig- 
inal down-spin electron trajectory has to be such that, 
as a result of the spin-orbit coupling, the electron spin 
reverts to the down position in respect to the direction of 
the field. The two consecutive spin flips imply the 2Ez 
energy dependence. 

At weak fields, when (2E Z / Esoi) 2 « T from 
Eq. (fT8|) . we obtain 



r,(0)rj(B) = 



E 9 



E. 



SOI 



(19) 



This is equivalent with 1 {B) — t < ^ 1 (0) = 
2 7 2 5 2 /(fcf,T 2 ), the result of Ref. [|. For 
(2Ez/E SO i) 2 - 1, Eq. (HID becomes: 



r.(0)rj(B) = -- 



2E, 



(20) 



1 - 



2E Z 
Esoi 



establishing the onset of the saturated behavior, that 
continues through the high field limit, when the dephas- 
ing rate is close to l/2r s (0) as observed experimentally 
0. In this regime, at a given intensity of the spin-orbit 
coupling Esoi > there is a maximum value of the longitu- 
dinal magnetic field determined by 2Ezmax — Esoi for 
which the modifications induced in the direction of one 
of the spins in the singlet Cooperon can be compensated 
by the skew-scattering produced by SOI, such that the 
singlet pairing is preserved. Increasing the intensity of 
the magnetic field beyond this limit, leads to the disap- 
pearance of the antilocalization correction as the singlet 
Cooperon configuration is not realized any more, exactly 
as observed experimentally [7j. The dephasing of the 
two trajectories remains constant since it is established 
solely by the spin-orbit coupling at a value equal to half 
the spin-relaxation rate along the z direction, (2t s (0)) _1 , 
2t s (0) being the amount of time in which the original 
down-spin flips and then realigns itself with the z axis. In 
the ideal 2D system modeled here, no surface effects have 
been considered. Consequently, this calculation does not 
reproduce the slow slope of the curve registered experi- 
mentally 0. 

We conclude, therefore, that a general theory of lo- 
calization effects in spin-polarized 2DEGs explains the 
saturated universal dependence on (2Ez/Esoi) 2 of the 
electron-dephasing time, matching the data of Meijer et 
al. 0. 
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